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The evolution of scalar linear perturbations is studied in gauge- invariant approach for 2-component 
models with nonrelativistic matter and minimally coupled scalar fields, the potentials of which were 
constructed for either constant dark energy equation of state (EoS) parameter w or its adiabatic 
sound speed equal to zero. The numerical solutions show that such fields are almost smoothed 
out on subhorizon scales. However they cause the scale dependent suppression of the nonrelativistic 
mater density perturbations and the decay of gravitational potential, which can be used for choice 
of the dark energy model. We discuss 2 types of the Lagrangian: classical and tachyonic ones. As 
our results show, the fields with w — const are almost indistinguishable, while for fields with = 
the difference of dark energy effective sound speeds c^, which is caused by the shape of Lagrangian, 
affects the evolution of perturbations significantly. We present also the transfer functions for both 
components. 
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I. INTRODUCTION 

Modern cosmological observations surely confirm that 
the expansion of the Universe is accelerated. Such phe- 
nomenon can be described using the models based on the 
Einstein equations with A-term [l|, 0] ■ Unfortunately, in 
this case several unsolved problems exist, e.g., the fine 
tuning and cosmic coincidence ones which are caused by 
the constant energy density corresponding to this funda- 
mental quantity Q. Alternatively one can assume the 
existence of new mysterious essence - dark energy. This 
repulsive component is regarded frequently as a scalar 
field, the Lagrangian of which can have either canonical 
kinetic term (classical field) 0, Q or non-canonical one 
(e.g., tachyon [6l|, k-essence [3]). 

It is often convenient to calculate the characteristics 
of the large-scale structure and CMB anisotropy in the 
phenomenological fluid framework, in which the dark en- 
ergy model is defined by its density, equation of state 
(EoS) parameter and effective sound speed. The dy- 
namical dark energy has been widely studied using both 
scalar field and fluid approaches on the background level 
as well as at the linear |8j, i9|] and non-linear [1^] stages of 
perturbations growth for either minimal or non-minimal 
coupling to other components. Moreover, the mod- 
els unifying both dark matter and dark energy have been 
proposed ^12;]. For quintessential dark energy the EoS pa- 
rameter w is assumed to be larger than —1, the models 
with w < —1 are called phantom ones [li|. The phan- 
tom divide crossing is not allowed for single fields with 
the simplest Lagrangians, therefore in such case the spe- 
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cific single-field models as well as multi-field ones should 
be discussed (e.g., [T^|). 

The purpose of this paper is to perform the compara- 
tive analysis of 4 possible scalar field models of dark en- 
ergy. This allows us to find the most sensitive to dark en- 
ergy type features and thus to propose the simple tests for 
identification of the source accelerating the expansion of 
the Universe. As the theoretically motivated scalar field 
potentials p^] frequently lead to the cosmological conse- 
quences contradicting the observations, we perform the 
reverse ingeneering, i.e., construct the potentials leading 
to the given dark energy properties (energy density, EoS 
parameter). We study the scalar fields with 2 shapes 
of Lagrangian: classical (Klein-Gordon) and tachyonic 
(Dirac-Born-Infeld) ones. For our research we have cho- 
sen the dark energy EoS parameter to be either constant 
or evolving in such way that the dark energy adiabatic 
sound speed equals 0. The reason for such choice is 
that in both cases the analytic solutions of the back- 
ground field equations exist. The latter parametrization 
is especially interesting since for homogeneous fields it 
produces naturally the behavior mimicing the dust mat- 
ter at early time and cosmological constant at very late 
epoch, so it can be regarded as possible candidate for so- 
called unified dark matter. All fields were studied assum- 
ing the same set of the best fitting cosmological param- 
eters obtained by WMAP and other projects. Since the 
background dynamics in models with both Lagrangians 
is the same, the main attention should be paid to the lin- 
ear stage of evolution of perturbations. On this level the 
degeneracy due to the Lagrangian type also exists, how- 
ever the wide class of models becomes principally dis- 
tinguishable. In this work we study the model of the 
Universe filled only with nonrelativistic matter and min- 
imally coupled scalar field acting as dark energy. We fo- 
cus on the infiuence of perturbed dark energy on matter 
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perturbations and emphasize the effect of the inclusion of 
dust matter on evolution of dark energy inhomogeneities. 
The latter effect is commonly ignored due to the unde- 
tectability of dark energy perturbations. Nevertheless it 
should be also investigated as it might be important for 
better understanding of the gravitational interaction of 
both dark components and thus for discovering of the 
nature of dark energy. We discuss also the gravitational 
instability in the single-component models with — 
in context of the unified dark matter. 

The paper is organized as follows. In Section |TT] we 
discuss briefly the background dynamics of our models. 
In Section IIIII the evolution equations for the gauge- 
invariant perturbation variables are written and adia- 
batic initial conditions are set. In Section llVI we discuss 
the numerical solutions of the perturbations equations 
for the single fields with — 0, 2-component model of 
the Universe and build the transfer functions for both 
components. We propose also the characteristics which 
can be used for determination of the dark energy type. 
The conclusions can be found in Section |Vl In Appendix 
we present the evolution equations for the field variables. 



II. COSMOLOGICAL BACKGROUND 

We consider the homogeneous and isotropic flat Uni- 
verse with metric of 4-space 



gi-jdx^dx^ — a^{ri){dr]^ — 5aiidx°'dx^), 



where the factor 0(77) is the scale factor, normalized to 
1 at the current epoch 770, rj is conformal time {cdt — 
a(r])dri). Here and below we put c = 1, so the time 
variable t = xq has the dimension of a length. The latin 
indices i, j, ... run from to 3 and the greek ones over 
the spatial part of the metric: a, /?, ...=1, 2, 3. 

If the Universe is filled with non-relativistic matter 
(cold dark matter and baryons) and minimally coupled 
dark energy, the dynamics of its expansion is completely 
described by the Einstein equations 



(m) 



(de) 



(1) 



where Rij is the Ricci tensor and Tj^p\ T^f*^^ - energy- 
momentum tensors of matter (tt^) and dark energy (de). 
If these components interact only gravitationally then 
each of them satisfy the differential energy-momentum 
conservation law separately: 



{ni,de) 
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(here and below ";" denotes the covariant derivative with 
respect to the coordinate with given index in the space 
with metric gij). For the perfect fluid with density p^m.de) 
and pressure P(m,de)i related by the equation of state 

P{m,de) = W(„i,de)P(m,de), it givCS 



P{m,de) — -'^-P{m,de)\^ + W(rn.de)) 



(3) 



(here and below a dot denotes the derivative with re- 
spect to the conformal time: "'"= d/drf). The mat- 
ter is considered to be non-relativistic, so Wm — and 
Pm ~ Pni a~'^ (here and below "0" denotes the current 
values) . 

We assume the dark energy to be a scalar field with 
either Klein-Gordon (classical) Lagrangian 



or Dirac-Born-Infeld (tachyonic) one 



(4) 



(5) 



where U{4>) and U{£,) are the field potentials defining the 
model, Xcias = ^;j^'V2 and Xtach = C;iC'V2 are kinetic 
terms. 

We postulate also that the background scalar fields are 
homogeneous (0(x, 77) = 0(f7), f(x,77) — ^(t;)), so their 
energy density and pressure depend only on time: 

Pclas = -^^'^ + U{(I)), Pclas = ^'/'^ " ^ {(f>) , (6) 



ptach = , ^^^^ , Ptach = ~U{0\n - ^- (7) 

The evolution equations for field variables ^(77) and ,^(77) 
are presented in Appendix \^ 

We specify the model of each field using 2 parameters: 
the EoS parameter Wde = Pdel Pde and the so-called adi- 
abatic speed of sound = pde/ Pde^ [E B [H, [i3> which 
satisfy the relation: 



Wde = 3aH{l + Wde)iWde " C^), 



(8) 



where H = d/a^ is the Hubble parameter (expansion 
rate) for any moment of conformal time 77. Generally, 
the equation of state is defined by the Lagrangian as 



L 



Wde 



(9) 



however, different shapes of L can lead to the same Wde 
- this is the well-known degeneracy on the background 
level. The scalar field evolution equations have the ana- 
lytical solutions for two cases: 

• W — const: = W, Pde{o) — Pde ^^"^ 

• = 0: w{a) = 7«oa^/(l + Wo ~ Wqu^), pde{a) = 
pf^[{\ + w,)a-^-w,]. 



^ The quantity doesn't play a role of the true sound speed of 
dark energy, it is rather a useful function, the form of which is 
similar to that of the adiabatic sound speed for barotropic EoS. 
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Here and below we omit index de denoting both - clas- 
sical and tachyonic - scalar fields for Wde- 

Note that the fields with = at a ^ 1 behave as 
the dust matter and at a — + oo they mimic the cosmo- 
logical constant. This fact suggests that they can also be 
studied as the possible candidates for the so-called uni- 
fied component, or unified dark matter, describing the 
early matter dominated stage and current A or dark en- 
ergy dominated one by introducing a single new fluid, 
e.g., Chaplygin gas [IJ]. 

In this case the expansion rate and fields evolve as 
follows: 



exist in each case [18|, |l9|. However, the physical con- 
sequences of both these solutions are the same, so from 
now we restrict ourselves only to the growing one. 

The analysis of dynamics of the Universe expansion 
for fields with w = const and = was presented in 
[Tsl [Toj . It doesn't depend on the scalar field Lagrangian 
and doesn't allow us to distinguish surely such models of 
scalar fields. So in order to propose the test for choice of 
the dark energy type which would be more adequate to 
observations we should study at least the linear stage of 
the evolution of scalar perturbations. 



H = Hoa "Sy/l + WQ- woa^,{10) 



/ y/l + ^;o(l " g^) - yi + wq 1 + yi + WQ ^^^^ 
\^yi + wo{l-a^) + y/l + wo 1- yi + Wo 



C(a) - Co = ± 



arctan 



' arctan 




For the 2-componcnt model the Hubble parameter 
H{a) behaves as: 



for w = const and 



3w 



H = Hoa 2 ^1 + QdeWo ~ ftdeWoa^ 



(13) 



(14) 



for cl = 0. 

The temporal dependences of scalar fields and their 
potentials are following: 



0(a) - 00 = ± 
U{a 



1 a'H{a') 
Pde{a) [1 - w{a)] 



for the classical Lagrangian and 

r°- da' y/l + w{a') 



^(a)-6 = ±^ 



a'H{a') 



U{a) = pde[a)\J-w{a) 



(15) 
(16) 

(17) 
(18) 



for the tachyonic one. 

In this paper we study the classical and tachyonic 
scalar fields with potentials constructed for w = const 
and = using the same set of the best fitting cos- 
mological parameters from {Vide = 0.722, w ~ — 
-0.972, a„ = 0.278, h = 0.697). For aU models there 
are 2 independent solutions for the field: the growing one 
(sign -|-) and the decaying one (sign -). Therefore 2 sym- 
metrical with respect to either (j) — 4>o or ^ — ^o potentials 



III. EVOLUTION OF SCALAR LINEAR 
PERTURBATIONS 

For derivation of the evolution equations for scalar lin- 
ear perturbations it is convenient to use the conformal- 
Newtonian gauge with space-time metric [20, 21] 



ds' 



a2(,y)[(l + 2^'(x,77))dr/2 
-{1 + 2^{:>i,Tj))Sapdx°'dx% 



(19) 



where ^(x, 77) and $(x, 77) are gauge-invariant metric per- 
turbations called Bardeen's potentials [13], which in the 
case of zero proper anisotropy of medium (as for the dust 
matter and scalar fields) have equal absolute values and 
opposite signs: \E'(x, 77) = — $(x, 77) [1^. In the linear 
perturbation theory it is convenient to perform the Four- 
rier transform of all spatially-dependent variables, so the 
equations will be written for the corresponding Fourier 
amplitudes of the metric (^'(fc,77)), matter density and 
velocity perturbations ((J^™-* (/c, 77), y(™)(A;,77)) as well 
as the scalar field perturbations {S'^'^'^\k,ri), V^'^''\k,r]), 
6(l){k,ri), 6^(k,ri)) (here k is wave number). These vari- 
ables are gauge-invariant [U, IH, IIJ] . The energy den- 
sity and velocity perturbations of dark energy, S^'^'^^ and 
yide) ^ are connected with the perturbation of field vari- 
ables S(f>, in the following way: 

S(^^-) = il + w)(^^-^ + ^^], (20) 



y(clas) ^ kScp 



(21) 



Sita.Hj^_l±^(SJ_^]^dll^^^ ^22) 



Ud^ 



y{tach) _ 



(23) 



Other non-vanishing gauge-invariant perturbations of 
scalar field are isotropic pressure perturbation 
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and intrinsic entropy 



(de 

r ^ 



^2 . 
Side 

W 



(26) 



The density perturbation of any component in the 
conformal-Newtonian gauge Ds = 5, which is gauge- 
invariant variable, is related to the other gauge-invariant 
variables of density perturbations D and Dg as: 



D = D„ 



, 3(1 + 
L», + 3(1 + w) 



w 

a V 
a k ' 



a k 



(27) 



where D^, D, Dg and V correspond to either m- or de- 
component. Here Dg is the density perturbation in the 
frame of vanishing curvature fluctuations and D corre- 
sponds to the rest frame, i.e., the frame in which the 
4-velocity is orthogonal to constant time hypersurface. 

For the scalar fields the intrinsic entropy perturbation 
is defined as [13, HI,!!!: 



(28) 



where is the dark energy effective (rest frame) sound 
speed: 



1 



dL 



Sp^'-f'> 



(29) 



V'^'^"^ + aH{l - ScDV^"^"^ - fc(l + 3cl)^ - 



elk 
l + w 



xDi'^"^ = 0, 



D. 



{de 



+ 3(c2 - w)aHDf''^ + (1 + w) 



k + -a^H^ (c- 



(33) 

2 



(34) 



Here and below p corresponds to the background density 
of each component. 

So, as we see, for each parametrization of EoS the 
shape of Lagrangian affects the evolution of perturba- 
tions only through the effective sound speed of dark en- 



ergy. In = consi-case c^^^ias) 



— 1' es{tach) — 



that 



in principle allows us to distinguish both Lagrangians, 
however the difference isn't large (for w close to —1 
- as it has been estimated on basis of the observable 
data 0) and suggests the similarity of solutions for 
both fields. In = 0-case c^^^j^^^ = 1 = const but 

e1{tach) ~ 7^ const, SO the behavior of perturba- 

tions in classical field and tachyon with the same cosmo- 
logical parameters should be really different. 

Therefore, in each case we have the system of 5 first- 
order ordinary differential equations for 5 unknown func- 
tions *(fc,a), Dl"'\k,a), V'^"'Hk,a), D'f''\k,a) and 
V'^'^'^\k, a). The following constraint equation is also sat- 
isfied: 



kH = AttGo' ( n„,D^^^ + pdeD^''- 



(35) 



For Lagrangians with canonical kinetic term it is always 
1, for tachyonic ones = —w. 

For further study it is convenient to use the evolu- 
tion equations for gauge-invariant density and velocity 
perturbation variables. The corresponding equations for 
perturbations of field variables can be found in Appendix 

El 

It should be noted that the principal degeneracy due 
to the Lagrangian shape exists also in the linear theory, 
i.e.jdifferent Lagrangians can lead to the same value of 
Cs However, the analysis of this stage of evolution of 
perturbations is still important since it removes partially 
the degeneracy existing on the background level. 



A. Evolution equations 

The linearised Einstein equations for gauge-invariant 
perturbations of metric, density and velocity are 



p„,V<^"'Hpde{l+w) 



■yide) 



= 0, 



D^ 



(m) 



(30) 
(31) 
(32) 



B. Initial conditions 

It is known that the observable large-scale structure 
has grown from the small adiabatic perturbations gener- 
ated in the early Universe. Since the density of the w — 
consi-fields is negligible at the early epoch (a <C 1) and 

= 0-fields mimic dust matter, all our models are ini- 
tially matter-dominated. In such case the growing mode 
of adiabatic perturbations corresponds to '^init = const. 
So, here we specify the adiabatic initial conditions for the 
growing mode of perturbations [23, 2^, : 



vide) _ ^ /T- 

^ imt — o LJ R T\ — y ^tmii 



^ init — " 



2 k 



3i?o VWI 



de 



3 init 

for w = const and 



yide) 



2 k 



init 



Jjide) 
9 init 



3 Hq VI + ^deWo 

2 k ^f.-^,, 



y/ Oiinit 1 



3 77 ^/YTV^^w^ 



^init 1 



(36) 
(37) 
(38) 
(39) 

(40) 
(41) 
(42) 
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FIG. 1: The Jeans scale for the Universe filled only with either 
classical (dotted line) or tachyonic (solid) field with — 0. 



9 init 



(43) 



for cl = 0. 

In all cases we have integrated numerically the sys- 
tems of 1-order equations for these initial conditions us- 
ing the publicly available code DVERK^. We assumed 
^init = —1, ainit = 10~^°. Such early initial time is 
chosen only in order to separate properly the growing 
adiabatic mode while our purpose is to analyze its be- 
havior after recombination epoch, when the Universe can 
be effectively described by the 2-component (dust matter 
plus dark energy) model. The evolution of perturbations 
is scale dependent, so we performed calculations for the 
range of k from 0.0001 to 0.1 Mpc~^. The reason for such 
choice is that at present epoch larger scales are super- 
horizon while smaller ones are affected by the non-linear 
effects. 



IV. RESULTS AND DISCUSSION 



A. Single-component system: perturbed scalar 
field 



First of all let us discuss the Universe filled only with 
a scalar fields. Such models with = 0, as the back- 
ground dynamics suggests, could probably be considered 
as simple candidates for the unified dark matter. There- 
fore we are going to analyze their clustering properties. 
For such purpose it is useful to obtain from the first-order 
equations using a second-order one for the 

density perturbation: 



3 

2a 



15 



2a 1 + wq 



H^a{l + wo — woa^) 



9 

^ 2^ 



D' 



WQa 



1 + Wq — Wga 



2 It was created by T.E. Hull, W.H.Enright, K.R. Jackson in 1976 
and is available at ,htt p : / / www .cs.t pronto. edu/NA/dverk.f.gz^ 



Woa 



3 
2^ 



D = 0. 



(44) 



1 -I- Wo — wofl 
The above equation could be written in the form: 

u" + uj^u = 0, (45) 

3 - 

where u stands for Da* (l + wq — wqq^) " . The coeffi- 
cient is generally a function of time and scale: 



H^a{l 
Woa'^ 



Wq 



woa^) 
27 



21 
16^ 



Wga 



1 + Wo — woa^ 16a^ \ 1 -I- — woa' 



However, in the neighborhood of any given point a it can 
be regarded as constant with respect to a, thus it depends 
only on the wave number k. The condition Lu'^{k) — 
specifies for each a the Jeans scale, under which the per- 
turbations at this a are stable (oscillatory) while above 
they are gravitationally unstable, i.e., can grow or decay: 



Hi r 



21 

Tea 



9wQa^ 



27woa'^ 



woa 



16 1 + Wq — waa^ 



(46) 



The Jeans scale for the values of a from 0.001 to 1 is 
shown in Fig. [T] For each field the values of k below 
the corresponding curve define the perturbations which 
can cluster for the scale factors close to a. As we see, 
the perturbations with k < 0.0001 Mpc~^ are unstable 
all the time up to the current epoch, however they are 
still super-horizon, so irrelevant for the choice of model 
best fitting to the data. The perturbations with subhori- 
zon scales could be unstable at early stages, but should 
oscillate at late times. 

In Fig. [5] the evolution of density, velocity and metric 
perturbations is shown. We see that generally the behav- 
ior of subhorizon modes is oscillatory. For classical La- 
grangian the gravitational potential decays very quickly, 
so that at a = 0.001 for k = 0.01 Mpc~^ (horizon entry 
at a « 0.0004) it is already almost while for k = 0.001 
Mpc~^ (horizon entry at a w 0.04) it is approximately 
—0.8 and starts to oscillate just before the horizon entry. 
For tachyonic Lagrangian the potential ^ at early stages 
is constant, then it begins to decay. Note that the per- 
turbations of tachyon with larger k change the character 
of the temporal dependence earlier (for 0.001 Mpc~^ at 
a w 0.09, for 0.01 Mpc"! at a w 0.03), but the amplitude 
of the first positive peak is almost the same (w 0.2). The 
gravitational potential for the super-horizon mode with 
k = 0.0001 Mpc""'^ decays monotonously, as expected. 

The rest frame pressure perturbations Sp'--'^-^^ — 
c^Sp^^^^ for cl > have the same sign as 5p^^^\ therefore 
the nature of the scalar field oscillations is similar to that 
of the acoustic ones in baryon-photon plasma. 

The presented in Fig. [5] evolution of density pertur- 
bations (related to as D = — fc^5'/(47rG'a^p)) confirms 
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FIG. 2: The evolution of the density (top) and velocity perturbations (medium) in a single-component model with = is 
presented for scales k = 0.0001 (dotted line), 0.001 (dashed), 0.01 Mpc~^ (solid). Bottom: the corresponding gravitational 
potentials. Classical field - left column, tachyon - right one. 



such scenario. We see that for the same scale the magni- 
tude of perturbation for the tachyon field is higher than 
for the classical one and that at late times all subhori- 
zon modes are smoothed out in oscillatory manner. The 
amplitudes of velocity perturbations oscillations at first 
grow slightly but then start to decay. 

The difference between the studied models is caused 
by the behavior of the effective sound speed. The per- 
turbations of the classical field always propagate with the 
speed of light (c^ = 1), so their evolution isn't close to 
that of dust matter at any stage of the growth. In con- 
trary, the perturbations of tachyon propagate with the 
speed variable in time: = —woa'^ / {1 + wq — woa^) . At 
early stages it is negligible, so such field mimics the dust 
matter not only at the background but at the linear stage 
of evolution of the perturbation too. Later, the effective 
sound speed grows and tends to the speed of light at in- 
finity (as well as the equation of state parameter tends 
to —1). The behavior of the field (and metric) perturba- 
tions begins then to differ from the dust matter one and 
soon becomes oscillatory with decreasing amplitude. 

Such features of the studied fields mean that they 
couldn't play the role of unified dark matter. 

The model with a single w = const-field corresponds 
to the quasi-de Sitter Universe. In this case the pertur- 



bation equation takes the form: 



9w 



2a 
~2~ 



2a 



. 3 
3w~- 



and has the analytical solution: 



D — a 1- Z 1 5+3u 

2 l + 3u 



„2i.2 



£1 = 



Cgk l + 3u 

-a 2 



l + 3w Ho 



(47) 



, (48) 



where Z^{x) = CiJ^{x) + C2Y^{x) and J^{x),Yi,{x) are 
Bessel functions of first and second kind, Ci , C2 - ar- 
bitrary constants. For a <C 1 the density perturbations 
behave as C'ia~2(i-™) -\-C2a~^~ - both modes decay for 
w < —1/3. In far future (a 00) the amplitude of D 
will decay as a 2 . Interestingly, the current oscillatory 
behavior of perturbations is clear also without analytical 
solution (|i5)) . Really, if we put equation (|T7)) in form of 
(|^ . the quantity 



UJ 



^0 



-,3w-l 



9u/_ 
16a2 



15w 

8^2 



21 
16a2 



is larger than up to present moment for all scales of 
interest. 
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FIG. 3: The evolution of variables D^'^^'^ {k = 0.01, 0.001 



-.{de) 



and 0.0001 Mpc"^ from top to bottom), D 
(corresponding scales from bottom to top) and V"''*'^) (k — 
0.001, 0.0001 and 0.01 Mpc"^ from top to bottom at a ^ 1). 
Classical field with w = const - dotted line, tachyonic one - 
solid. The curves for both fields are very close and practically 
indistinguishable. 



matter clusters while dark energy tends to homogeneity. 

For the dark energy with w = const the simple conclu- 
sion, that the behavior of scalar linear perturbations in 
model with the tachyonic field should be similar to that in 
the model with the corresponding classical field, is valid, 
as the numerical analysis has shown |28(|. In plots pre- 
sented in Fig. [3]it is hard to distinguish the curves for the 
classical field and tachyon, so in this case we practically 
can't choose the Lagrangian preferred by observations. 
For scale k — 0.0001 Mpc^^ the behavior of dark energy 
density perturbations is different: Dg'^'^^ remains almost 
constant while D^'*'^) and o'f''^^ grow, but this depends 
only on gauge choice for super-horizon modes. 

In the = 0-case both fields are almost smoothed 
out on subhorizon scales at present and future epochs 
while on superhorizon ones they do not grow significantly. 
However, the difference of the effective sound speed be- 
havior has sufficient imprint in the evolution of inhomo- 
geneities, so, as it can be seen in Fig. |31 the perturbations 
of the classical field at early stages grow slowly but decay 
after the horizon entry while the perturbations of tachyon 
at first grow significantly and then begin to oscillate. The 
oscillation amplitudes decrease all the time for the den- 
sity perturbations or increase at early time while decrease 
at the current epoch for the velocity ones. It should be 
noted that small oscillations along the averaged solutions 
for and D'^'^"^ are also present in classical field, how- 
ever their amplitudes are highly subdominant comparing 
to the corresponding mean quantities. Interestingly, for 
perturbed tachyonic field at scale k — 0.001 Mpc~^ Dg 
remains positive. This suggests that the small-scale per- 
turbations of tachyon oscillate along the mean curves too, 
but in this case the averaged values are much smaller (al- 
most negligible) comparing to the amplitudes of oscilla- 
tions. 

The infiuence of dust matter component on dark en- 
ergy perturbations is described by the equation: 



B. 2-component system: dust matter and dark 
energy 

Now we return to the minimally coupled 2-component 
model. The general conclusion is that in the studied case 



where 



(49) 



J 



A 
B 

S 



^ , fc2(3^_3^2_i)^127rGaVn(3w;-3c2-2)\ ^ 
ow H — — — TT-^ I ati. 



k^ + UnGa^pm 
. fc2(3w - 3cl + UTrGa^prniSw - 3cl - 2) 



(50) 



P-,12.GaV. 3.(aH)--9a^i/^(l + .)(.- a 

47rG'aVde(l + w) -l- 12nGa'^ p,ncl, (51) 
D^-^ + 3aH fs. -3cl-2- " + UnGa^i^u, - 3c^ - 2)^ 



fc2 + l2TTGa^p„ 



xinGa'^Pmi^ + w). 



(52) 
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FIG. 4: The evolution of density perturbations D^'^'^\ D^^'^\ Di'''^\ velocity ones V^'^'^^ and gravitational potential with scales 
k = 0.0001 (dotted hne), 0.001 (dashed) and 0.01 Mpc~^ (solid) for classical (left) and tachyonic (right) fields with — 0. 



For p„j = it takes the form of ([33|). It can be easily 
seen that for the same values of cosmological parameters 
the non-zero density of nonrelativistic matter changes the 
coefficient B and thus the character of evolution of the 
dark energy perturbations more drastically for fields with 
larger effective sound speed. 

In the bottom panel of Fig. |4]the evolution of —4' is 
shown for models with = for the scales of pertur- 
bations k = 0.0001, 0.001 and 0.01 Mpc'^ It can be 
easily seen that in this case the strong model and scale 



dependences of ^' exist. 

Generally, at present epoch the growth of matter den- 
sity perturbations (described by the growth function g) 
is suppressed and - unlike ACDM-case - such suppres- 
sion is scale dependent. The growth function is defined 
as follows: 

.g^i?(™)a,„,/i?(:::la. (53) 

It is convenient to compute this quantity using the solu- 
tions of (|30p - (IMl) and taking into account (PT)) . In the 
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FIG. 5; The evolution of g for linear perturbations with scales 
k = 0.0001, 0.001, 0.01 and 0.1 Mpc"^ (from top to bottom 
in each panel) in the models with non-relativistic matter and 
dark energy with w = const (top) or = (bottom). Clas- 
sical field - dotted line, tachyon - solid line. Both axes are 
logarithmic. 



top panel of Fig. [5] the temporal evolution of g is shown 
for the fields with w = const. We see that all curves prac- 
tically overlap. It means that the scale dependences are 
very weak and confirms the absence of substantial differ- 
ence between classical and tachyonic dark energy. The 
temporal dependences of g for the fields with = are 
shown for different scales in the bottom panel of Fig. [51 
The suppression caused by the classical field is stronger 
than the corresponding suppression caused by tachyon. 
The behavior of the growth function can be used for 
choice of the adequate model of dark energy because 
it is connected with the dust matter transfer function 
and thus with the observable matter density fluctuations 
power spectrum. 

The faster decay of ^E" and caused by the classical 
field, is clear. From the equation 



^- + 3-* 

a 



2^ 
a 



a k 



(54) 



we see that for models with the same EoS parameter and 
adiabatic sound speed of dark energy the influence of its 
perturbations on the gravitational potential and thus on 
dust matter inhomogeneities growth increases with the 
values of the dark energy effective sound speed. As for 
classical field always equals 1 and for tachyon it ini- 
tially is close to while at the current epoch it equals 
—Wo J which is smaller than 1, the perturbations of the 
former field affect the matter and metric perturbations 
more. The scale dependences of ^ and g, caused by 
the dark energy perturbations and therefore absent in 
ACDM-cosmology, are also expected to be stronger for 
classical Lagrangian. 



It should be noted that the oscillatory behavior of per- 
turbations of the tachyonic field leads to small oscillations 
along the averaged solution for 4", but their amplitudes 
are subdominant comparing to the decaying mean depen- 
dence. 

The dependences of the gravitational potential and 
growth function on scale of perturbations for the present 
epoch are shown in Fig. [H] The presented plot supports 
the conclusion that the scale dependence in w = const- 
case is very weak and almost the same for both types of 
the Lagrangian (as mentioned above, a bit higher sup- 
pression of g caused by classical field is due to the ef- 
fective sound speed, which in this case is slightly larger 
than for tachyon). At small scales \1/ and g are prac- 
tically scale-independent - here the behavior is close 
to that in the "quasi-ACDM" model, i.e., model with 
w = const 7^ —1 and without dark energy perturba- 
tions ^28^. For c^ — the scale dependence for classical 
field is significantly stronger than for tachyonic one (in 
agreement with the result deduced from temporal depen- 
dences). The oscillations seen in the scale dependence of 
^' for model with the tachyonic Lagrangian simply reflect 
the fact that small oscillations along the mean temporal 
dependence of gravitational potential come to a = 1 for 
different k in different phases, so they lie above or below 
the averaged curve. 

In far future (a ^ cxd) all perturbations, which in 
past have entered the particle horizon, will exit the event 
one. Both fields with = tend to the cosmological 
term mimicry, therefore the perturbations will be entirely 
smoothed out. In w = const-case the situation is more 
complicated. When the inhomogeneities in dark energy 
are neglected, all perturbations will decay to 0. However, 
taking them into account leads to the freeze of gravita- 
tional potential after the event horizon exit, as it follows 
from the numerical analysis and asymptotic solutions. 
The dark energy perturbation variables D'^'^^'^ and T/^'^*^' 

decay but d'^^^^ and d'^^^'^ tend to the constant values. 
It should be mentioned that the latter variables decay 
in sub-horizon regime but after the event horizon leav- 
ing they start to grow in order to reach the frozen val- 
ues. Nevertheless, the described features are irrelevant 
for choice of the best fitting dark energy model because 
they correspond to the future behavior of super-horizon 
and thus non-observable inhomogeneities. 



C. Transfer functions 

The observations of the large-scale structure of the 
Universe provide us with the power spectrum of mat- 
ter density perturbations P(fc, a), which is related to pri- 
mordial (post-inflationary) one via the transfer function 
r(fc, a) as follows: 

P{k,a) = Ask'''T^{k,a), 

where Ag is the normalization constant, - scalar spec- 
tral index. 
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FIG. 6: The scale dependences of gravitational potential (top) and growth function (bottom) at the present epoch for dark 
energy with w = const (left) and = (right). Classical Lagrangian - dotted line, tachyonic one - solid. 



The transfer function can be built for each component 
of the Universe, and for dark energy too [26]. In this pa- 
per we are interested in scale dependences of the transfer 
functions 



Tm{k,a) 



(55) 



and 



D(^^){km.n,a) D(<^^)(fc,a™„) ■ ^'^^ 

Here amin is a very small scale factor, for which all scales 
of interest are still superhorizon. In Fig. [7] these de- 
pendences are presented for the present epoch (a — 1). 
It can be seen that in the w ~ const-case the scale de- 
pendence of is weak and both curves overlap. It is 
an illustration of the fact that the influence of perturba- 
tions of such fields on matter ones is almost negligible. 
For comparison we show the cold dark matter (CDM) 
transfer function computed for w — const using CAMB^ 
[29| . Such calculation takes into account not only the 
gravitational coupling of the dark energy to pressureless 
matter but all physical processes governing the evolution 
of its perturbations. We see that the influence of dark 
energy inhomogeneities on CDM ones is almost negligible 
comparing to the effect caused by other components and 
their interactions. 

In = 0-case the transfer functions decay more 
quickly than in w = const-one. This fact suggests that 



^ http://camb.info 



the influence of such fields on the computed by CAMB 
CDM transfer function (corresponding to the real Uni- 
verse) could become comparable to the other processes 
determining it. It should be noted that in agreement with 
previous conclusions the perturbed classical field affects 
the matter transfer function more than tachyonic one. 

In the bottom panel the scale dependence of Tj^ is 
shown for the present epoch. As we can see, compar- 
ing with the dust matter transfer functions the dark en- 
ergy ones are reduced by approximately 9 orders for both 
w — const-fields, 10 orders for the classical field with 
= and 5 orders for = 0-tachyon at small scales 
(0.1 Mpc~^). The scale dependence of Tj^ is monotonous 
for w = const-fields with both Lagrangians and = 0- 
classical field (the latter is smaller approximately by an 
order than the former ones). For tachyonic field the scale 
dependence of transfer function oscillates with decreas- 
ing amplitude. These oscillations are simply produced 
by different phases, which the density perturbations of 
different scales have at a = 1 . It is interesting to see that 
their period (especially for large k) is almost constant 
and equals approximately 0.00045 Mpc~^. 



V. CONCLUSION 

The evolution of the scalar linear perturbations is stud- 
ied for the 2-component (dust matter and minimally cou- 
pled dark energy) cosmological models. The dark energy 
component is assumed to be either classical or tachyonic 
scalar field with the potentials constructed for EoS pa- 
rameters satisfying the condition either w — const or 
cl = 0. 

The analysis of gravitational instability of the single 
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FIG. 7: The transfer functions for dust matter (top) and dark 
energy (bottom) at a = 1. The computed by CAMB cold dark 
matter transfer function for fields with w = const is shown 
for comparison. 



speed of dark energy. The discussed scale dependence 
provides us with simple tests for choice of model which 
describes the observational data in the best way. The 
tests can be based on measurements of the matter den- 
sity perturbations power spectrum and late ISW-effect. 
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APPENDIX A: EVOLUTION OF SCALAR FIELD 
VARIABLES AND THEIR PERTURBATIONS 



field models shows that the behavior of subhorizon per- 
turbations is oscillatory in all studied cases. Thus the 
fields with = 0, although on the background level they 
mimic dust matter at early epoch and the cosmological 
constant in far future, can't serve as unified dark matter 
because they don't cluster on scales smaller than particle 
horizon. 

In all 2-coniponent models the dark energy is almost 
smoothed out on subhorizon scales. The perturbations of 
fields with w = const grow or decay without visible os- 
cillations and their evolution is similar for both types of 
Lagrangian (Fig. [S]). The behavior of the ~ 0-tachyon 
density perturbations remains oscillatory (unlike that of 
the corresponding classical field, see Fig. [4]). Such differ- 
ence is caused by behavior of the effective sound speed of 
dark energy. In w = const-case the quantity is large 
and almost the same for both types of Lagrangian at all 
stages of evolution of perturbations as soon as the ob- 
servational data prefer the values of the EoS parameter 
close to —1. In = 0-case the effective sound speed of 
classical field is always equal to 1, while of tachyon 
varies in time as —w, so it is relatively small before the 
beginning of dark energy dominated epoch. It was found 
that the inhomogeneities of scalar fields with larger ef- 
fective sound speed are more affected hy pm ^ This 
effect, though undetectable, could be interesting from the 
purely theoretical point of view as it might be useful for 
better understanding of the dark energy nature. The 
same applies to the dark energy transfer functions. 

The caused by perturbed dark energy scale dependence 
of the nonrelativistic matter growth and transfer func- 
tions as well as of the gravitational potential (Fig. [5][7l) is 
small and almost the same for both fields with w = const. 
In = 0-case it is large and different for classical field 
and tachyon. Such behavior is due to the effective sound 



Evolution equations for the background scalar fields 
0(77), ^(77) and their linear perturbations 5(j){k^ rj), S^{k, rj) 
can be obtained either from the Lagrange-Euler equa- 
tions or from the differential energy-momentum conser- 
vation law Tr 



{de 



0:i 



= 0. 



In the discussed cases they are following: 



2aH^ + a'^^ = 0, 



<j) + 2aHS(f) 
-4«'0 



+ a- 



6(p + 2a' 



.dU 



for fields with Klein-Gordon (classical) Lagrangian and 




5^ 



for fields with Dirac-Born-Infeld (tachyonic) one. 
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